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ABSTRACT 



We present the results of a search for molecular gas emission from a star- 
forming galaxy at z = 4.9. The galaxy benefits from magnification of 22 ±5 x due 
to strong gravitational lensing by the foreground cluster MS1358+62. We target 
the CO (5-4) emission at a known position and redshift from existing Hubble Space 
Telescope/ ACS imaging and Gemini/NIFS [Oll]3727 imaging spectroscopy, and 
obtain a tentative detection at the 4.3 ex level with a flux of 0.104±0.024 Jykms -1 . 
From the CO line luminosity and assuming a CO-to-H 2 conversion factor a = 2, 
we derive a gas mass M gas ~ llo.6 x 1O 9 M . Combined with the existing data, 
we derive a gas fraction M gas / (M gas + M*) = 0. 59^06- The faint line flux of 
this galaxy highlights the difficulty of observing molecular gas in representative 
galaxies at this epoch, and suggests that routine detections of similar galaxies in 
the absence of gravitational lensing will remain challenging even with ALMA in 
full science operations. 

Subject headings: galaxies: high-redshift — galaxies: star formation — gravita- 
tional lensing: strong 



1. Introduction 



Extensive studies of optically-selected star-forming galaxies seen at the epoch of peak 
cosmic star-formation density (z ~ 2) have reve aled star-formation rates (SFRs) of 10 



lOOMfpyr 1 and stellar masses of ~10 10 11 M (e.g. iForster Schreiber et al.ll2009l ; lLaw et al. 



20091 : iTacconi et al.ll2010l : iDaddi et al.ll2010bl ). If these SFRs have been continously main- 



tained, then these galaxies must have undergone their first major epoch of stellar mass 
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assembly 1-2 Gyr earlier, at z ~ 5, when the bulk of the star- forming population was ~ 5x 
less massive (jMcLure et al.ll2009l : Istark et al.lboioh . 



The star formation within these galaxies is fuelled by reservoirs of predominantly cold 
molecular hydrogen, H 2 . Since the H 2 is not directly detectable, CO emission at millimetre 
wavelengths has been employed to trace the cold molecular gas. However, exploring gas 
properties of galaxies beyond z ~ 3, is challenging; not only does the apparent surface 
brightness reduce as (l + z)™ 4 , but the galaxies themselves also appear systematically smaller 
and intrinsically fainter making detections of their molecular gas emission difficult. 

To date, studies of molecular line em ission at z > 3 ha ve been limited to extreme 
populat ions such as gas-rich quasars (e.g. IWalter et al.l 120091) and submillimetre galaxies 



(SMGs; ICoppin et al.ll2010l ; ICarilli et allBOld iRiechers et al.ll2011af ). Detecting typical star- 
forming galaxies at z > 3 has proven difficult b ecause the CO line l uminosities are usu ally 
below the sensitivity limits of current facilities (jStanway et al.ll2008l ; iDavies et al.ll2010l ). 



However, it is still possible to study less active high-z galaxies which are strongly grav- 
itationally lensed by massive galaxy clusters. The magnification provided by gravitational 
lensing has e nabled detections of molecular gas in a number of star-forming galaxies up 



to z ~ 3 (e.g. iBaker et al.l 12004 ; ICoppin et al. 



201ll ). and more recently z ~ 4 — 6 (e.g. ICox et al. 



2007 



Riechers et al. 



2011 



20 



Combes et al. 



lb; Danielson et al. 



2012|). The physical 



properties of the interstellar medium appear similar to those in local ULIRGs, with high gas 
fractions, high densities and intense UV radiation fields. 



Franx et al.l (119971 ) reported the detection of a multiply-imaged z = 4.9 galaxy which is 



gravitationally lensed by the massive galaxy cluster MS 1358+62. Correcting for lensing, this 
galaxy (hereafter MS1358-arc) appears to be representative of the star-forming population at 
this epoch (its lensing-corrected apparent magnitude is Iar =< 24.9, margin al ly brighter than 



the ch aracteristic luminosity of Iab = 25.3 at z ~ 5; lOuchi et al.l I2004J ) . ISwinbank et al. 
( 120091 . hereafter S09) carried out a detailed study of one image of the MS1358-arc system, 
using optical and infrared imaging combined with integral field spectroscopy, revealing a 
rotating system across 2 kpc in projection, with star formation occuring in five bright clumps. 

In this Letter, we report observations with the Plateau de Bure Interferometer to search 
for CO (5-4) emission in MS1358-arc. Throughout, we adopt a ACDM cosmology with Hq = 
70kms" 1 Mpc" 1 , Q A = 0.7 and Q m = 0.3. 
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Observations and Data Reduction 



The MS1358+62 system is illustrated in Figured! with the positions of the two images 
of MS135 8-arc marked (Table U). We estimate t he amplificat i on of each image using the lens 
model of iRichard et al.l ( 120081 ) with lenstool ( IKneibl Il993l ; I Julio et al.l 120071 ). with errors 
accounting for the magnification gradient across the image. The combined magnification 
factor for Image 1 and 2 is /x = 22 ± 5. 

To search for the CO (5-4) emission, we observed MS1358-arc with the IRAM Plateau 
de Bure Interferometer (PdBI; iGuilloteau et al.lll992l ). using six antennae in D (compact) 
configuration. The pointing centre was ct^otxb^ooo = 13 : 59 : 48.7, +62 : 30 : 48.34, and 
the frequency was tuned to the CO (5-4) transition (rest frame 576.2679 GHz) redshif ted t o 
97.185 GHz (based on the [Oll]-derived systemic redshift of z = 4.9296 ± 0.0002 fromH). 
The observations were made between 2010 June 4 and 2010 June 6 with total on-source time 
of 10 hours, using WIDEX with a resolution of 2.5 MHz. The star MWC349 was observed 
as the primary flux calibrator, with the quasars 1749+096, 0923+392 and 2145+067 as 
secondary flux calibrators. Receiver bandpass calibration was performed against 0923+392 
and 3C454.3, and 1435+638 and 1418+546 were used for phase and amplitude calibration. 



The data were reduced with the gildas software package ( IGuilloteau fc Lucad 120001 ) 
and resampled to a velocity resolution of 38.6 km s -1 . The synthesised beam is a Gaussian 
ellipse with a FWHM size of 5.6 x 4.3" at a position angle of 111°. 

The PdBI observations were centered on the two brightest images of the z ~ 5 galaxy. 
The images are both within 6" of the pointing centre, whereas the primary beam of PdBI 
at 3mm is 50"; we have therefore not corrected for primary beam attenuation, which is 
negligible. First, we construct a channel map from a 90km s -1 window around the CO(5- 
4) emission line centered at the systemic redshift, and show as an inset in Figure [TJ with 
contours at ±lcr intervals. The positions of the images within the PdBI cube are determined 
by aligning the cube with the HST image and taking the centre of the corresponding pixels. 

We also extract spectra from these positions and show these in Figure [2j Individually, 
the two images are clearly extremely faint; by comparing the A% 2 of a Gaussian profile to 
that of a continuum-only fit, measuring the noise over a 1200kms _1 channel, the significances 
of the CO (5-4) emission lines in Images 1 and 2 are 3.3cr and 2.7cr respectively. However, 
since they are separated by more than one beam width, we can improve the signal to noise by 
coadding the spectra at the two positions using uniform weights (Figure |2]); we then obtain 
a la rms channel sensitivity of 0.4 mJy per 38.6kms _1 channel. Again comparing the \ 2 °f 
a Gaussian profile fit to that of a continuum-only fit, we use the A% 2 to derive a significance 
level of 4.3cr for the coadded CO (5-4) emission line. 
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Fig. 1. — Hubble Space Telescope (EST) / ACS image of the core of MS1358+62. The squares 
mark the extent of the pixels in the PdBI cube used to extract the spectra shown in Figure |2j 
A Plateau de Bure channel map constructed from a 90km s _1 window around the systemic 
velocity of the CO (5-4) line is shown as an inset. The PdBI synthesised beam size is indicated 
by the dotted ellipse in the inset. Contours are shown in ±ler intervals, with negative 
contours indicated by dashed lines. The origin of the HST image coordinate system is 
(«2ooo^2ooo) = (13 h 59 m 51.4',+62 o 30'52.5"). 
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Table 1: Positions of the three images of MS1358-arc within the Plateau de Bure field of view 
and the mean linear magnification factor (fj,). The position of the third image is included 
for completeness, but is excluded from our analysis as it lies outside the primary beam. 
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Fig. 2. — The CO (5-4) spectra extracted from the two images of MS1358-arc and the coadd 
of these spectra. The velocity scale is centered on the systemic redshift z = 4.9296, and 
the expected position of the CO (5-4) emission line is indicated by a black dotted line. The 
spectra are binned into 38.6kms _1 channels. The best-fit Gaussian profile to the coadded 
spectrum is overlaid as a guide, scaled to the relative amplifications o f the two images, and 
the velocities of the two brightest star-forming clumps measured by |S09| are indicated by 
vertical lines. Each image shows positive flux at the systemic redshift of MS1358-arc, but 
individually they are only detected at ~ 3a significance. 
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Whilst the emission line centroid is extremely well-matched to the best-fit redshift of the 
nebular emission, the line is weak, highlighting the difficulty of these observations even with 
long integrations. We perform a number of tests to validate the robustness of the detection. 
First, we check how often a detection is made when coadding two random spectra from 
the original (non-coadded) cube. We generate 1000 combinations of two randomly-selected 
pixels on the map (excluding pixels which lie within one beam of MS1358-arc), then fit a 
Gaussian profile to the resulting coadded spectrum and compare this to a continuum-only 
fit. The probability of finding an emission line of equal significance to the target from two 
randomly-selected positions with a velocity centroid lying within ±150kms~ 1 of the redshift 
of MS1358-arc and a line width in the range a = 30 — 250 km s -1 is 0.05%. 

The CO (5-4) emission line we measure is centered at z — 4.9297 ± 0.0001 with a line 
flux of 0.104 ±0.024 Jy km s" 1 and Full Width at Zero Intensity of FWZI = 150 ±20 km s" 1 , 



comparable to the velocity gradient measured across the source by IS09I . Error bars are 
obtained by fitting a Gaussian profile to the line and perturbing it in intensity, velocity and 
width to obtain a A x 2 = 1 error surface, using a Monte Carlo routine with 10 5 realisations 
centered on the best fit. 

We also construct a coadded channel map from the PdBI cube by extracting regions cen- 
tered on the two images and co-adding them. The coadded cube is then spatially convolved 
with the PdBI beam, and spectrally convolved with a Gaussian with FWHM = 100 km s -1 . 
We estimate the noise at each spatial position using the off-line spectrum at that position, 
and divide the signal in the on-line slice of the convolved datacube by this noise map to 
construct a signal-to-noise map, shown in Figure[3l We find that 99.95% of the pixels have 
lower signal-to-noise than the target. This is equivalent to a detection significance of 3.5 a. 

In the following analysis, we treat the line flux as 0.104 ± 0.024 Jykms" 1 , but discuss 
the implications of treating it as an upper limit in Section HI 



3. Results, Analysis and Discussion 

3.1. CO Dynamics, Luminosity and Molecular Gas Mass 

Using the combined spectrum of Images 1 an d 2, w e first estimate the line luminosity 
and gas mass. Following ISolomon fc Vanden Boutl (120051 ). the line luminosity is -^co(5-4) = 
(3.6 ±0.8) x 10 9 Kkms _1 pc . Adopting a magnification factor of ji = 22 ± 5, this indicates 
an intrinsic ^co(5-4) = (1-6 ± 0.5) x 10 8 Kkms^ 1 pc 2 . 

To convert this CO (5-4) luminosity to a total molecular gas mass, we must first esti- 
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Fig. 3. — Upper: The coadded PdBI signal-to-noise channel map, constructed by extracting 
and coadding regions around each of the two images of MS1358-arc, convolved spatially with 
the beam and in the spectral direction with a Gaussian profile with FWHM= 100km s _1 , and 
then dividing each pixel by its noise. The combined image of the two targets is located in 
the centre of the map. The solid lines are intensity contours in ±lcr intervals, with negative 
contours indicated by dashed lines. The beam size is indicated by the black ellipse line in 
the lower right-hand corner. Lower: Cumulative histogram of the signal-to-noise in the map, 
with the central pixel marked as a red dashed line. We find that just 0.05% of pixels have 
higher signal-to-noise than the emission seen from the source, suggesting a 3.5<r detection. 



-8 - 



mate the corresponding CO (1-0) line luminosity. We therefore exploit recent observations 
of z = 2 — 4 galaxies where the ratio L' CQ , 5 $\/L' c o n ) has been measur ed directly, e.g. 

^CO(5-4)/£'cO(l-0) 



0. 35 ± 0.02 for SMM J2135 flDanielson et all 1201 If ) or 0.32 ± 0.05 



for a sample of SMGs (IBothwell et al.ll2012l ). We adopt the latter value and applying it 

(5.0 ± 1.7) x 10 8 Kkms 1 pc 2 . If we instead as- 



to MS1358-arc, we estimate L[ 



sumed L 



CO(1-0) 



CO(1-0) 



from eit her the Antanna e, 0.2 (IZhu et all 120031 ): NGC253, 0.7 



(IBavet et al.ll2004f ): or M82, 0.7 f lWeiB et all 120051 ) this estimate would vary by a factor 
of 2x in either direction and so we adopt this as the likely systematic uncertainty in our 
estimate. 



To derive the molecular gas mass A^f gas uum ^00(1-0) ' vvc arouiuc -^^gas — ^^cofi— 0) 



from L' co ^_ j, we assume M gas = aL' c 
where a is the coefficient relating the CO to H 2 +He gas mass. We can derive an absolute 
lower limit for a by assuming that the molecula r gas i s enriched to solar metallicity and is 
optically thin to CO radiation. From llvison et al.l (120101 ). we have a > 0.7, which is consistent 
with the value a = 0.8 applie d to local (Ultra) Luminous InfraRed Galaxies ((U)LIRGs 



Solomon fc Vanden Boutll2005l ). However, recent studies of high-z star-formi ng; galaxies have 
suggested that a _~ 2 m ay be a more appropriate conversion factor (e.g. iDanielson et al. 



20 111 ; llvison et al.l 1201 lh . Adopting this value, we find M gas = ll 06 X 10 9 M Q , where we 



conservatively estimate a factor of 2x uncertainty in a, comparable to the uncertainty in 



J CO(5-4) 



/L'c 



V CO(1-0)- 

Combining our gas mass estimate with the stellar mass M* 
by IS09I . we obtain a total baryonic mass of M baryon = M* + M gaf 



7 ± 2 x 10 8 M derived 



1.7+J;° x 1O 9 M . This 
3 ± 1 x 



is consistent with the dynamical mass within 2kpc derived by IS09I of Md yn 
M , if the inclination is % > 45°. 



10 9 csc 2 (i 



The gas fraction / gaf! = M gas / (M gas + M*) is typically < 10% in local larg e spiral 



galaxies (lYoung fc Scovilld Il99lh . or ~ 33% in local ULIRGs (ISolomon et al.lll997l ). Other 
studies of molecular gas in hi gh-redshift galaxies are beginning to fin d a trend for higher gas 



fractions at higher redshifts ( ITacconi et al.ll2010l; iGeach et al 



expectations from hydrodynamical simulations ( ICrain et al.ll2009f ). For MS1358-arc, we find 



20 1 lh . which is in line with 



/gas = 0.59io;Qg, with a systematic uncertainty of ~ ±0.20 due to a and L 



CO(5-4) 



/L'c 



v CO(l-0)' 

In Figure HI we co mpare the gas fraction of MS1358-arc to samples at lower redshift (see 
also lGeach et al.ll201ll ). The sample size is currently too small to draw any firm conclusions, 
and there are of course significant selection effects involved. We therefore caution that there 
are limits to the extent of any physical interpretation of this result, although it clearly 
motivates a uniformly-selected survey of the gas properties of high-redshift star-forming 
galaxies. 
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We can also use the dynamics of the [On] emission to test whether the molecular gas is 
colocated with the star formation. In Figure [21 we overplot the redshifts of the two brightest 
star-forming knots, which are located at ±150 km s" 1 from the dynamical centre. Hence, if 
the CO gas traced the [Oil] emission we would expect a FWZI of ~ 300km s _1 , higher than 
the observed line width. This may indicate that the CO emission is associated with the 
dynamical centre of the galaxy rather than being concent rate d in the star-forming regions 
or associated with the outflowing Lya and UV-ISM lines ( 1S091 ). 



3.2. Gas depletion timescales 



S09l derived a star formation rate (SFR) of 42 ± 8 M© yr" 1 using the iKennicuttl (119981 ) 
conversion from [On] luminosity to SFR. W e note that thi s cony ersion is consistent with 
the empirical correction to this calibration of iGilbank et al.l ( 120101 ). which accounts for the 
metallicity dependence of the [Oll]-derived SFR. 

If MS1358-arc has sustained this SFR continuously, then the time taken to build up the 
current stellar mass of 7 ± 2 x 1O 8 M is just Tbuiid ~ 17Myr (consistent with the age and 
star formation history used to derive the stellar mass). This suggests we may be seeing this 
galaxy in its first epoch of star formation. If this SFR is maintained at a constant level, 
the gas supply of 1 x 10 9 M will be exhausted in ~ 24 Myr, which would place the galaxy 
~ 40% of the way through its starburst. The star formation 'lifetime' of this galaxy would 
thus be very short compared to that inferred indirectly for the LBG popu lation at z = 4 — 6, 
which have implied starburst lifetimes of < 500 Myr (IStark et al.l 120091 ). but is consistent 
with the lifetime of th e starburst-triggered 'LBG phase' predicted by semi-analytic models 
jGonzalez et al.lboilh . 



Finally, we note that the ratio between far-infrared luminosity and molecular gas mass, 
Lfir/Mh 2 , gives the star formation efficiency (SFE), relating the ongoing star formation to 
the available molecular gas. 



Knudsen et al.l (120081 ) obtained an upper limit on the 850 /im flux of MS1358-arc of 
Ssso < 4.8 m Jy which, accoun ting for the lensing magnification and using the Ssso — -^fir 



calibration of iNeri et al.l (120031 ) . suggests a far-infrared (FIR) luminosity Lfir < 3.5 x 10 11 L Q 
(we note that if we assu me that the FIR fl ux is re-emitted light from dust-obscured star 
formation and apply the IKennicuttl (119981 ) relation to our [Oil] -derived SFR, we obtain 
-^fir = (2.3 ± 0.9) x 1O 11 L , which is consistent with this upper limit). We thus obtain 
SFE < 24OL0MQ 1 . For comparison, the median SFE in local ULIRGS is 49 ± 6LQM" 1 
when corrected to a = 2, although this rises to 172±23L Mq 1 if a = 0.8, which is the factor 
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Fig. 4. — Evolution of the molecular gas fraction with redshift in star-f orming galax- 
ies. We compare MS1358-arc to the z = samp l e of iLerov et al.l (120081 ). the median 
gas fractions of the high-;? LIRGs fro m iDaddi et al.l (]2010al ) and the median of each red- 
shift b in of the iTacconi et al.l ( 120101 ) (U)L IRG sample. We also use the lensed ULIRG 



from Kneib et a 



Bothwell et al.l ( 120121 ) 



(2005 ), lensed LBGs from iRiechers et al.l (2010al ) and high- z SMGs from 



Riechers et al.l ( 12 01 Obi ) and Walter et al.l (120121 ). all converted to 



a = 2 for consistency. The lines are the predictions from semi-analytic models for galaxies 
with SFR > 1, 10, 100 M Q yr _1 (Lacey et al. in prep). From the limited data available, there 
is apparent evolution in the molecular gas fraction with redshift, broadly consistent with the 
semi-analytic models. 
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commonly applied to ULIRGs flSanders et al.lll99lh . while SMGs give SFE ~ 55^L M" 1 



for a = 2 (IBothwell et al.ll2012f ). Our upper limit is consistent with these values. 



4. Conclusions 

We have used the IRAM Plateau de Bure Interferometer to search for the CO (5-4) 
emission in a z = 4.9296 galaxy lensed by the foreground cluster MS1358+62. We were 
able to observe two images of the galaxy simultaneously, with a total magnification factor 
of 22 ± 5x. We measure a line flux of 0.104 ± 0.024 Jy km s _1 at the position of the galaxy, 
yielding a detection at 3.5-4.3cr. 

The molecular gas shows a relatively narrow velocity range around the systemic redshift, 
unlike the [Oil] emission, which mainly arises from two clumps in the galaxy at ±150 km s -1 . 
This may suggest that the gas is more centrally concentrated than the star formation. 

We derive a total gas mass M gas = ll 06 x 1O 9 M , suggesting that this galaxy has 
a gas fraction / gas = 0.59^^06' w hh a systematic uncertainty of ~ ±0.20 due to a and 
-^co(5-4)/-^co(i-o)' wri ich is similar to the most gas-rich galaxies at z ~ 2. This could imply 
that gas fractions do not continue to rise significantly beyond z ~ 2, though a larger sample 
is clearly needed to draw any conclusions. 

Finally, given the tentative nature of the detection, we consider the implications of 
treating the measured flux as a 4cr upper limit. In this case, the resulting gas mass M gas < 
1 x 10 9 M Q would be lower than expected given the stellar and dynamical masses derived by 



S09_|, and the gas fraction would be / gas < 0.6. This would also imply a gas depletion timescale 



of < 24Myr, placing the galaxy more than 40% of the way through a short starburst. 

Our observations highlight the difficulty of measuring gas properties of 'representative' 
star- forming galaxies at z ~ 5, an era when many of today's massive galaxies may be 
undergoing their first major episode of star formation. Probing their basic properties - 
their stellar and gas content and relation to star formation - will provide important physical 
quantities which galaxy formation models must reproduce. 

The authors thank Malcolm Bremer, Seb Oliver and Isaac Roseboom for help and 
useful discussions. RCL acknowledges a studentship from STFC AMS acknowledges an 
STFC Advanced Fellowship, and RGB, IRS and ACE acknowledge support from STFC. 
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